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1. Introduction 

Here I return to the previously considered [1] random wandering ( "Brownian motion" ) of 
a particle interacting with ideal gas. Investigation of this system, with the aim to obtain 
complete actual quantitative statistical characteristics of the Brownian motion, still 
remains an intriguing mathematical problem for the future. The principal preliminary 
analysis of this problem, as well as of microscopic Brownian motion in fluids in general, 
was undertaken in [2, 3, 4, 5, 6, 7] and even much earlier in [8, 9] (with continuation in 
[11, 10]), and in respect to various non-fluid systems in [12, 13, 14, 15] and [9]. Basing 
on exact complete sets of evolution equations for many-particle probability distribution 
functions, - the Bogolyubov-Born-Green-Kirkwood-Yvon (BBGKY) equations [16], - 
it was argued that long-scale statistics of the Brownian motion qualitatively differs 
from the Gaussian one, thus rejecting the "law of large umbers". This means that 
diffusivity and mobility of a molecular Brownian particle (BP) have no definite values 
but instead undero low-frequency fluctuations with some 1/f-type spectrum. Although 
that are natural thermal fluctuations, they were lost by conventional "Boltzmannian" 
approaches to fluid kinetics, since these approaches neglect many-particle statistical 
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correlations. But the 1/f fluctuations appear as soon as one takes into account all the 
infinite chain of many-particle statistical correlations. Paradoxically, the latter originate 
from unpredictability of dynamical cause- and- consequence correlations between inter- 
particle collisions, as it was explained first in [12, 13] and [8] and later in [9, 11, 2, 4, 5, 7]. 
Most mathematically rigorous results in this field are presented in [3] and [4] while most 
developed approximate solutions to the BBGKY equations in [8, 9, 11, 10]. 

Nevertheless, there is a lack of formally simple and at the same time irrefutably 
convincing proofs of the foregoing statements. Here I suggest such a proof. It covers 
not only thermodynamically equilibrium Brownian motion but also (and first of all) 
non-equilibrium one under a constant external force applied to the BP. 



2. The equations and the problem 

The BBGKY hierarchy describing BP in ideal gas were written out in various 
representations and discussed in [1, 3, 4]. Therefore, please, see firstly Sections 2 and 
3 in [1]. Then now we can start from the evolution equation for generating functional 
[3, 4, 16] of the cumulant correlation functions: 
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where P = MV and p = mv are Momentum = Mass*Velocity of the BP and gas atoms, 
respectively, R is BP's position (and its total path if it starts from the origin R = 0), 
p represents relative distance of atoms from BP, x = {p, p} , f . . . = f . . . dx , v is 
mean gas density, $(p) is BP-atom (repulsive) interaction potential, 

G m {p) = (2nTmy 3 / 2 exp(-p 2 /2Tm) , 

E(r) ee exp[-$(r)/T] , (2) 
E'(r) =dE(r)/dr = — $'(r) E(r)/T , 

and, at last, if)(x) is arbitrary bounded function, and / is the external force. In fact, 
presence of this force and corresponding operator / • d/dP (entering any equation of 
the BBGKY hierarchy) is the only difference of Eq.l from equations of [1]. 
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(3) 



Notice also that 

CG M (P) = , (5) 

hence, Gm{P) is equilibrium stationary solution to Eq.l at / = . 

As before in [1], let us assume that at initial time moment t — the BP was 
disposed at definitely known point R = while the gas was in equilibrium (that is the 
force / is switched on just at t — ). Then initial condition to Eq.l again is 

V{t = 0} = G M (P)5(R) , (6) 

Solving Eq.l with this condition, we can find the probability density distribution of BP's 
variables, W = W(t, R, P \ v, f) , - which eventually is most interesting for us, - from 

W = V{ip = 0} 

At ip 7^ , the functional V = V{t, R,P,ip\ v, /} represents full statistics of BP's 
motion in a "cloud" of its correlations with surrounding gas atoms. 

Importantly, as it was underlined in [2, 3, 4], the virial relations [1, 2, 3, 4, 5, 6] 
remain valid in presence of the external force. Thus at any t, R, P and / we can write 

V{a + i> | u} = V{^/(1 + a) | (1 + a) v} (7) 

with a = const being arbitrary constant from interval — 1 < a < oo . One can easy 
derive this "generating virial relation" directly from Eq.l and condition (6) by exploiting 
clearness of the form (3). This relation, (7), will play a crucial role below. 

For further it is convenient to go from V to corresponding characteristic functional 
V defined by the Fourier transform 

V = J J exp(ik-R + i£-V)VdPdR , 

and to the BP's characteristic function 

W = J f exp(ik ■ R + i£ -V) W dPdR = 
= = 0} = 

= {exp[ik.R(t)+it.V(t)] \uj) , 
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where we introduced angle brackets (. . . \u, f) to denote statistical averaging at given 
v and / under mentioned initial conditions. 

The angle brackets' designations help to represent statistical contents of V and V 
in very transparent form. Namely, in accordance with the V 's definition [1] % , 



Here 



_ _ ( exp [ik ■ R(t) + i£ ■ V(t) + / In (1 + ji{x)) v(t, x)]\u, f) 
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is microscopic gas density in the \i -space, and 
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thus expressing the Poissonian statistics of ideal gas in equilibrium. 
The Eq.l and the initial condition to it transform into 



V{t = 0} = exp[-T£ 2 /2M] 



Si/j(x) 



+ vg(x) 



(8) 



(9) 

(10) 
(11) 



The problem is an adequate analysis of this equation from viewpoint of the characteristic 
function W . 



3. Cumulant representation 

Next, it is reasonable to introduce cumulants and their generating functional: 



V = exp S , 

S = 1^ a \ h \ S ab {t^\vJ} 

a,b = 

Evidently, 

(i) any of the functionals S a f,{t,ifj \ u, f} , at a + b > , represents all mutual 
irreducible correlations of (a + b + n) -th order between a samples of BP's path R(t) , 
b samples of BP's velocity V(t) and n samples, - with n — 1, 2, . . . , - of instant gas 
state v(t, x) . At n — , S a b{t, \ v, /} is (a + b) -th order mutual cumulant of R(t) 

t 

See also [2, 3, 4] and, for details, the basic Bogolyubov's definition [16]. 
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and V(t) irrespective to the instant gas state. Taking into account that under our 
initial condition R(t) = f Q V{t')dt' , and denoting purely irreducible correlations (i.e. 
cumulants) by double angle brackets, we can write 

S ab {t, 0|i/,/} = ((R a (t)V\t)\u,f)) = (12) 

f ((V(t 1 )...V(t a )V b (t)\u,f))dt 1 ...dt a ; 
o 

(ii) because of the relation (7) any of the functionals S ab {t, ip\v,f} possesses 
similar property: 

S ab {t, a + i> \u, /} = S ah {t, ^/(l + a)\(l + a) u, /} (13) 

with arbitrary — 1 < a = const < oo ; 

(iii) in terms of the functional S = S{t,k,£, ip\v,f} Eq.9 takes the form of 
nonlinear equation with quadratic nonlinearity, 




dt~ k d£ +f M 



+ 



J[i + 1,( x )]l(^ + 



S{t = 0} = -T£ 2 /2M 



d_ 



% 1 

M 



5S 



+ vg(x) 



(14) 



(15) 



Of course, eventually we are most interested in the generating function of BP's 
cumulants themselves, S{t,k,£, | is, /} , which gives diffusivity, mobility and other 
statistical characteristics of BP's motion. But anyway we should start from realizing 
some principal properties of the whole set of cumulants contained in S{t, k,£, ip\ u, /} . 



4. Time evolution and spatial extension of BP-gas cross-correlations 

Although probability distribution of the BP's path R(t) constantly evolves with time, 
distributions of BP's velocity V(t) and gas state v(t,x) can be expected (at least, 
under proper interaction potential) to tend to a stationary limit. Then 

V{t, k = , ^, i> 1 1/, /} ->■ V stat {£, ip\v,f} = cxp S stat {£, ip\v,f} , 

where the limit function obeys Eq.9 at k = and zero time derivative: 

= / ~ V-tat + C' Vstat (16) 

Let us discuss this equation. 

The structure of the operator C ' , as well as of C , obviously allows solution of 
Eq.16 be extended from any "good" ip(x) (e.g. vanishing at x — > oo ) to ip(x) + a , - 
with a being a constant (— l<o"<oo),-by equating 

Sstatti, a + V I v, /} = S stat {i, V/(l + o) I (1 + <t)v, /} (17) 
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Combining this virial relation with its analogue (what directly following from (7) or 
(13)) for any finite t , we have 

S s tat{Z, <? + i> k, /} - S{t, 0, £, a + v k, /} = (18) 



where 



^'(s) = ^ , 1/'= (l + a)u 

1 + (T 



The equality (18), in spite of its seeming triviality, says, together with (17), about 
several important things as follow. 

Statement 1 . 

All stationary cross-correlation cumulamts 

Sfi'(x 1 ...x n \v,f) = ((V b v(x 1 )...v(x n )\uJ)) 

are different from zero and (absolutely) integrable functions of Xj . Thus all they vanish 
when any of pj or pj goes to infinity. In other words, all the BP-gas correlations 
are localized in the p -space near BP, and corresponding characteristic "volume of 
correlation" [2, 3, 4] is finite. 

In essence, this property of the stationary cumulants, - as well as analogous property 
[1, 2, 3, 4] of the arbitrary non- stationary ones, 

S ahn (t, Xl ...x n \u,f) = ((R a (t)V b (t)v(t,x 1 )...v(t,x n )\v,f)) , 

- is necessary boundary condition for a correct construction of solutions to Eq.16 and 
the evolution equations. 

Statement 2. 

Choosing ip(x) in Eq.18 to be a "good" function, well localized in the p -space, 
we can conclude that local values of non- stationary cumulants and their integrals over 
some of Xj tend to their stationary limits with one and the same speed. 

In other terms, none of the cumulant functions Sobn{t, X\ . . . x n \u, f) can contain 
such a component, c(t, x\ . . . x n \ v, f) , that 

->■ (19) 

In particular, for example, Sobiit, x\u, f) can not have a part, c(t,x), behaving 
somehow like 

g(x)e(ut-\p\) 



c(t, X\ . 


■■X n W,f) 


J c(t,xi . . 


.X n \u,f) dXj 



c(t,x) oc 



->■ , J c(t,x) ->■ const ^ , (20) 



with 6{-) being the Heaviside function and u some characteristic velocity. Such 
(contributions to) correlations can be named "phantom correlations", for they 
simultaneously "thaw" at infinity and stay significant, leaving nonphysical "invisibly 
small BP's correlations with infinitely far points of gas". 

Hence, Eq.18 (and its parent, Eq7) rejects "phantom correlations", thus producing 
even more strong restrictions on possible evolution of the correlations than their spatial 
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integrability, and ensuring finiteness of the "volume of correlation" during all the 
evolution. 

Statement 3. 

The two above statements naturally must be extended from S bn(t, x 1 . . . x n \u, f) 
to general cumulants S a b n (t, x i ■ ■ ■ x n \u, f) including BP's path R(t) . 

Indeed, firstly, the exclusion of the "phantom correlations" (possessing (19)) merely 
determines more rigid (and physically meaningful) boundary conditions at infinity 
than the spatial integrability in itself. Secondly, although some of the cumulants 
S a b n (t,xi . . .x n \v, f) at a > definitely have no time limits, all they concern the 
same asymptotically stationary state of the system as cumulants with a = . Hence, 
they are determined by the same boundary conditions. 

Notice that mathematically exclusion of the "phantom correlations" means 
commutativity of limits t — > oo and i[>(x) — > const , while physically finiteness of a 
number of gas atoma actually involved into BP-gas correlations. 

Let us apply this to the simplest cumulants Sio n (t, x\ . . . x n \u, f) and their 
generating functional Si (t, tp \ u, f) , and substantiate one more principal 

Statement 4 ■ 

Correlations between the total BP's path R(t) and instant gas state v(t, x) 
generally (at / ^ ) are growing with time under the same law as the path itself. 

Clearly, far enough in the stationary state, at / ^ , mean value of BP's path must 

grow proportionally to time, ((R(t))) = ((V(oo))) t + const , while cross-correlation 

cuulants, {{R(t) v(t,Xi) . . .v(t,x n ))) , either tend to constants or grow not faster than 

((R(t)}) (otherwise the Statement 4 is even more than true). Therefore, firstly, we can 

introduce limit 

c (i)/ , i f \ - y Sio(t,ip\v,f) 
SioWWJ) = hm 

t ^,oo t 

Secondly, the above Statement 3 about the boundary conditions at infinity allows us to 
extend virial relations (13) to the limit functional, writing 

S$>(* + i/,\v,f) = S$>U,/(l + a)\(l + *)u,f) (21) 

and 

* + (22) 

= 5!'VK/)- S ' o(i '^ K/) ^o 

(withe above defined primed variables). 

The limit virial relation (21) unambiguously claims that the limit functional 
Siq (ip \u, f) is actually depending on ip(x) , since otherwise (21) would claim 
independence of the stationary mean BP's drift velocity lim ((R(t)))/t = ((V(oo))) 
on the gas density v (which is obviously unacceptable). Hence, all the limits 
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lim ((R(t) v (t, Xi) . . . v (t, x n )))/t are non-zero finite quantities. This is the same as 
the Statement 4. 

Simultaneously we have came to 

Statement 5. 

The time-averaged drift velocity, R(t)/t , is actually random quantity regardless of 
duration of the averaging. 

Indeed, otherwise all the limits lim ((R(t) v(t, x±) . . . v(t, x n )))/t (n > 0) would 
be equal to zero. 

Next, we will consider fluctuations of the drift velocity. 

5. Real long-time asymptotics versus conventional model asymptotics, and 
BP's mobility fluctuations 

Excluding somehow from the evolution equations (1) or (9) or (14) the gas related field 
variable ip(x) we would obtain a closed but time-nonlocal equation for BP's variables or 
cumulants.. Such operation can be named "exclusion of thermostat" or "derivation of a 
kinetic equation for BP". In fact, the "exclusion of thermostat" never was realized in a 
honestly correct way. Instead, it always was based on a priory neglect of inter-particle 
statistical correlations (at least, three-particle and higher-order ones). The result of this 
assumption can be e.g. the Boltzmann-Lorentz equation [17]. 

Unfortunately or fortunately § , an unprejudiced investigation shows || that in 
reality all the many-particle correlations always are significant. To feel once again 
what does it mean, we have to compare a long-time asymptotics of BP's cumulants (12) 
dictated by Eq.14 with the asymptotics following from model kinetic equations, e.g. the 
Boltzmann-Lorentz equation. 

Firstly, let us recall the standard 

5.1. Model asymptotics 

Kinetic equation for BP in homogeneous media looks like 



where C mo d e i is kinetic operator, - generally integral one and certainly satisfying 
£-modei Gm{P) = (analogue of our equality (5)), - for instance, the Boltzmann-Lorentz 
operator (i.e. linearized Boltzmann operator). Corresponding equation for 





I think "fortunately" since otherwise the world would be too primitive and boring! 



See e.g. [8, 4, 11, 5 



7] and especially the Krylov's prophetical book [18]. 
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with the same initial condition (15). 

Consider asymptotic behavior of the cumulants generating function, S(t, k, £ \ v, f) , 
at t/r — > oo , where r is characteristic momentum relaxation time . The asymptotic 
has the well known and clear form, 

S(t,k,t\v,f) = SW(k\v,f)t + S®(k,t\u,f) + ... , (23) 
where the dots replace remaining terms what are decaying to zero, and 
c(D h dSi0) _, f * r f dS(0) a. 9 

This linear asymptotics means, obviously, that second- and higher-order irreducible self- 
correlations of BP's velocity are fast decaying functions (decaying exponentially or at 
least in an integrable fashion) of time differences, so that all cumulants (12) with b > 
tend to finite constants, while at b = correspondingly grow proportionally to time. 

In other words, - borrowed from the probability theory [19, 20], - R(t) 
asymptotically behaves as a random process with independent increments and thus 
has asymptotically Gaussian probability distribution. 

Now, let us return to our exact equations and perceive that 

5.2. Exact equations forbid the model asymptotics 

In ideal gas, regardless of its density, certainly there are no collective excitations and 
thus no hydrodynamical correlations. Therefore, for the first look, we can expect 
that the Boltzmann-Lorentz equation gives qualitatively correct description of BP's 
motion, moreover, even quantitatively correct description of its asymptotical statistical 
properties. In other words, we would like to expect that at t/r — > oo , similar to (23), 

S = t + 5 (0) + . . . (24) 

At that, can be more formally introduced by = lim S/t . Again, such the 
asymptotics would mean that all irreducible self-correlations of BP's velocity, as well 
as its mutual correlations with gas state, v(t,x) , are fast enough decaying (integrable) 
functions of time differences. 

This expectation, however, immediately meets serious objections and contradic- 
tions. We will consider them in a few steps. 

Notice that the asymptotics (24) can not be literally analogous to (23) since factor 

must be essentially dependent on ip . Indeed, from the exact relation (7) or (13) it 
follows, - in view of arbitrariness of t and in view of results of the previou Section 4, - 
that must obey similar relation, 
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S^{k,^,a + ^\u,f} = S^{k,^,^/(l + a)\(l + a)u,f} , (25) 

with all potential arguments being written out. This relation allows be independent 
on £, like S w in (23), but requires its actual dependence on ip . Otherwise (excluding 
tp from the list of its arguments), we inevitably would come to conclusion that is 
also completely independent on v . This, in turn, would imply that mean value, variance 
and higher cumulants of R(t) , and hence BP's diffusivity and mobility, are completely 
independent on gas density ! 

Since the latter is certainly wrong statement (i.e. such S can not bring solution to 
Eq.14), we have to accept expression (24) in the form 



S = SW{k,tl>\u,f}t + SW{k,Z,rl>\v,f} + ... (26) 
with being a non-trivial functional of ip(x) , so that 

5ip(x) ' Sip(xi) 5vjj(x2) 
(the sequence is infinite, because v(t,x) by its sense is non-Gaussian random field). 
Thus we came to 

Contradiction 1 . 

From (26) and (27) it follows, for large enough t , that 

({R(t)v(t,x)\u,f)) cx t , (28) 
((R a (t)v(t,x 1 )...v(t,x n )\is,f)) oc t , 

where all (omitted) coefficients on the right-hand sides generally (at / ^ ) are 
essentially finite, that is non-zero, non-negligible and non-vanishing with time. In 
particular, in respect to the first row of (28), 

r 1 {{R(t)v(t,x)\isJ)) -+ c(x\uj)t , (29) 

where the (restored) coefficient c(x\v, f) is finite in the above sense, so that 

c(x\v,f)^0 (30) 

Hence, mutual irreducible correlations (cumulants) of BP's velocity and gas state, 

{{V(t 1 )v(t,x)\is,f)) , 

(( V(h) . . . V(t a ) v(t, Xl ) . . . v(t, x n ) \v,f)) , 

are not fast decaying (integrable) functions of t — tj . This in hard contradiction with 
what was assumed as a ground for (24) and (26) ! 

Clearly, attempt to include £ into list of arguments of can not improve the 
situation. The matter is that just emphasized non-integrability of of mutual, or cross, 
correlations of velocity and gas state, t.g. ((V(ti)v(t,x) \u, /)), necessarily implies 
non-integrability of velocity's self-correlations, e.g. (( V(ti) V(t) \ u, f )) , and, as the 
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consequence, violation of the assumed asymptotics (24). This will be seen soon in 
Remark 2. 

Remark 1 . 

Considering infinitesimal form of relation (25), one can easy obtain 

((R(t)v(t,x)\u,f)) = u^-((R(t)\u,f)) , (31) 

which is example of particular virial relations [1]. Combining it with (29), we have 

c(x\u,f) = lim u — ((V drift (t)\u,f)) , (32) 

where V dri f t (t) = R(i)/t is time-averaged velocity, or "drift velocity", of BP. 
Remark 2 . 

Expressions (29)- (32) together do prompt that the "drift velocity" Vd r ift{t) is 
essentially random quantity even at arbitrary long duration of time averaging. The 
word "essentially" underlines that magnitude of fluctuations of Vd r ift(t) is comparable 
with its mean (ensemble average) value. This statement follows already from the (exact!) 
relation (31) if supplemented with reasonings expounded in [2, 3, 4] % . 

We can come to the same statement merely if combine the widely known general 
inequality (in essence, the Cauchy-Buniakowski inequality) 

((AB)} 2 < ({A 2 }} ((B 2 )} 

with (29). Let A = R(t) and B = J v(t, x) x(x) , with x( x ) — In [1 + ^(x)] being 
some suitably fixed function (see (8)). Then 



X (x){(R(t)v(t,x)\isJ)) 



< 



(33) 



< (( 



X(x)v(t,x) 



u,f))((R\t)\u,f)) 



Applying (29), we have 



«[/*(*) «(*,*)] I"./)) 

In view of (30) and (32), evidently, x( x ) here always can be chosen such that 
the coefficient before t 2 is non-zero. Hence, the variance of the drift velocity, 
a/ (( R 2 (t) | v, f )) / 1 , is bounded from below by a finite quantity independent on the 
observation (averaging) time ! 

Thus we came to 

1 

A volume in the p -space actually contributing to left side of (31) has the order of l/v (at least, 
when BP's mass is comparable with atom's mass). 
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Contradiction 2. 

The inequality (34) says that the second row of (28) is in in contravention of the 
first row. Thus inequality (34) is principally incompatible with the assumed hypothetical 
asymptotics (24) what produces both the rows simultaneously. 

Eventually, we arrive to 

Conclusion. 

The model asymptotics (24) is forbidden by the exact Eq.l (or Eq.9 or Eq.14) as 
it was claimed in the title of the present Section. . 

This means that true statistics of the Brownian motion qualitatively differs from 
statistics implied by conventional model kinetic equations. The difference manifests 
that real BP has no certainly predictable drift velocity and mobility (or, in other words, 
its mobility possesses slow "quasi-static" fluctuations). 

6. Magnitude of the mobility fluctuations 

The same conclusion would appear if we tried to search for solution of Eq.14 in the form 
(26). But it would require much more tremendous consideration than the above one 
based on the virial relation (7). Nevertheless, of course, any quantitative calculations 
of statistical characteristics of the Brownian motion are impossible without direct 
investigation of the evolution equation (14) (or (1) or (9)). In particular, calculation of 
the important cross-correlation cumulant function c(x\i>, f) , - introduced in (29), - as a 
functional of the BP-atom interaction potential (or corresponding "scattering matrix"). 
Since, however, methods for solving of such functional PDE as (1) or (9) or (14) still 
are not developed + , we are forced to confine ourselves by rough estimates. 

Therefore, returning to inequality (34), let us choose x( x ) — <f>(p) ■ Then (34) 
changes to 

(( R 2 (t) I *,/»> [MpK/MpW 2 e > (35) 



where 



and 



(([Jd>(p)u(t,p)dpY\u,f)) 
J (pWJ) = J c(x\v,f)dp 
<(t,p) = J v{x)dp = 5(p-(rj(t)-R(t))) 



is random microscopic gas density in the configurational space. According to (32), 

c'(p\u,f)dp = lim u—((V dnft (t)\u,f)) (36) 

+ 

Although some formal approaches were suggested in [1] and approximate ones in [8, 11]. 
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To estimate the coefficient before t 2 in (35), firstly, let us assume that the external 
force is sufficiently small, e.g. in the clear sense that ( = \f\ VqTq/T <C 1 with 
Vq = yjT/M , while the observation time is sufficiently large, e.g. in the sense of 
t ^> tq/( 2 (which means that drift component of the BP's path R(t) is mach greater 
than its diffusive component* ). Then denominator in (35) can be estimated as 

2 



J <t>(p)Ht, P ) 



dp 



« v J E{p)(t) 2 {p)dp w i/ J <p 2 (p)dp , 

while Vdriftit) and c'(p\u,f) in (36) and R(t) in (35) represented as 
VdHftit) = pit) f , 

R(t) = p{t) ft , 
c'(pWJ) = vp\p\v)f , 

where p(t) plays the role of random "small field" mobility. After that (35) reads 

{(p 2 {t)\v)) = (p 2 (t)\v)-(m\") 2 > 



> 



v[J P'(pW)(f>(p)dpf 
fo 2 (p)dp 



while (36) turns into 

p'{p\v) dp 



d_ 

dv 



(mi*) 



(37) 



(38) 



(in respect to first-order cumulants double and single brackets are equivalent). 

Secondly, let us choose in (37) such 0(p) what maximizes the right-hand expression. 
The maximization yields 



(p 2 (t)\v)-(p(t)\v) 2 > » I P'\p\v)dp 



(39) 



The rest of estimate is less formal. Let £1 = be a characteristic finite volume 
(space region) in the p -space, such that it produces main contributions to the integrals 
in (38) and (39). Then, obviously, 

[J n p'dp] 2 [jp'dp] 2 



p' 2 dp > I p ri dp > 



.'2 



(40) 



Besides, notice that from physical point of view the only natural measure for the 
characteristic volume is the specific volume 1/V . Therefore we have rights to write 
il(v) ~ 1/V . Adding these reasonings to (38) and (39), we find 

1 2 



(p 2 (t)\v)-(p(t)\v) 2 > 



(41) 



Formulas and pictures for "molecular Brownian motion" under such regime were presented and 
discussed e.g. in [21, 22]. 
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A rigorous version of this estimate can be formulated as follows. For any < a < 1 
let Q(a, v) be the minimum of all volumes (space regions) satisfying 



/ p' dp — I p' dp 
Jn J 

Then one can verify that 



< a 



p' dp 



> (I -a) 



p' dp 



p' dp 



Applying this inequality in place if the last step in (40), we have 



p dp > 



(l-a) 2 [jp'dpf 



Q(a, v) 

Finally this yields, instead of (41), 

(p 2 (t)\v)-(p(t)\") 2 > 



1 2 



:i-^) 2 



max 



(42) 



(43) 



This correction, however, does not cancel the estimate (41) Jj . 

We see that magnitude of the mobility fluctuations generally can be expected on 
order of its mean value. 

Remark 3. 

Here, to keep logics, we should answer two questions as follow. 

Just made estimates exploited the finiteness (expressed by (30)) of the function 
c(x\v, f) which in turn had appeared as a part of the hypothetical asymptotics (26). 
But the latter was logically rejected! Then why we can use one of its consequences? 
And why one can not imagine this function to be tending to zero with time while its 
support, Q , in the p -space growing to infinity in such way that its integral in (32) 
stays constant? Under such scenario the estimate (43) would become insignificant. 

In fact, the answers already were done in Section 4. 

Firstly, asymptotics (26), being wrong as the whole, at the same time is true 
in respect to the first-order terms in ik -expansion of the S{t, k, £, /} , more 
precisely, in respect to Sio{t, /} . Secondly, the above imaginary scenario involves 
nonphysical "phantom correlations" and therefore is forbidden by the virial relations. 
All this forms sufficient ground for our estimates. 



7. Discussion and resume 



To conclude, let us point out most principal aspects of our above consideration. 

(i) The exact evolution equation (1) for generating functional of many-particle 
cumulant correlation functions of the system "molecular Brownian particle (BP) in gas" 
[1, 4] produces exact "virial relations" [1, 4] connecting various statistical characteristics 

tt 

About heuristic physical reasons for the estimate Ct(v) ~ 1/v sec [2, 3, 4]. 
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of the system, in particular, in the above considered non-equilibrium steady state driven 
by an external force applied to the BP. 

(ii) The virial relations, in their turn, imply obligatory requirement to all irreducible 
(cumulant) correlations between BP's velocity and total path, from one hand, and 
current microscopic gas state, from the another hand, to be always (even in the non- 
equilibrium steady state) integrable functions of distances between BP and gas atoms. 
Besides, all these cross-correlations must be located near BP, never running away to 
infinity ff . 

(iii) As the consequence of the locality of BP-gas correlations, magnitude of BP's 
path and gas state cross-correlation and variance of the path are on the same order of 
value as the path itself, irrespective to the time of evolution and path observation. 

This means that BP's mobility has no certainly predictable value but instead 
undergoes slow quasi-static fluctuations whose magnitude is comparable with its 
ensemble-average value. 

Thus, the exact evolution equation, in opposite to various model kinetic equations, 
predicts existence of BP's mobility 1/f noise. 

(iv) The integrability of BP-gas correlations says that the latter envelope only a 
finite number of gas atoms (on order of unit) in BP's vicinity. Hence, neither BP nor gas 
remember a history of their interaction (conserved in states of far running away atoms). 
Consequently, the system has no possibilities to "control and regulate" a number of 
BP-atom collisions and thus BP's mobility. 

This is just those reason of the mobility 1/f fluctuations what for the first was 
guessed in [12, 13, 8]. 

Notice, besides, that the inequality (43) can be treated as "uncertainty 
relation" between mean square of of the mobility fluctuations and the "correlation 
volume". At that, treating the latter as a measure of of the system's memory about 
its past, we come to statement not once pronounced in my cited works: the shorter is 
system's memory, the greater is its 1/f noise. 

(v) Importantly, all our above consideration can be easy generalized from BP 
in ideal gas to BP in arbitrary fluid, if using the results of [3, 4] (then, for instance, 
generalization of the estimate (43) will differ from (43) by additional multiplier at its 
right-hand side, Tdu/dV , where V is gas pressure). 

But, besides the logical analysis of the evolution equations ("Bogolyubov equations" 
[4]) and principal estimates of their solutions, we are interested also in potentially exact 
regular methods for analytical solving of these equations (in addition to approximate 

tt 

In this respect it is interesting to notice that truncations of the BBGKY hierarchy leading to model 
kinetic equations (and thus losing the mobility fluctuations) simultaneously born non-local running 
away correlations which involve unbounded number of atoms. 

Such correlations may look like the "phantom correlation" (20). In correct theory, it would transform 
into something like g(x) exp (—\p\/X)/p 2 T . This subject will be considered separately. 
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methods from [8, 11, 10, 21] and formal "boson representation" from [1]), in order to 
calculate spectra and probability distributions of the mobility fluctuations. 
This may be very difficult but intriguing adventure. 
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